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Abstract We analyze a series of full-Sun observations, which was performed with the SoHO/SUMER instrument between 
March and October 1996. Some parameters (radiance, shift and width) of the S vi 93.3 nm, S vi 94.4 nm, and Ly s line profiles 
were computed on board. Radiances and line-of- sight velocities in a large central region of the Sun are studied statistically: 
distributions of solar structures, field Fourier spectra and structure functions are obtained. The structures have distributions with 
power-law tails, the Fourier spectra of the radiance fields also display power laws, and the normalized structure functions of 
the radiance and velocity fields increase at small scales. These results support the idea of the existence of small scales, created 
by turbulence, and of intermittency of the observed fields. These properties may provide insight into the processes needed 
for heating the transition region, or, if confirmed in the corona, the corona itself. The difficulties encountered in this analysis, 
especially for the velocity data, underline the needs for sensitive ultraviolet imaging spectrometers. 
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1. Introduction 

The mechanism heating the solar corona to millions of degrees 
kelvin remains an open problem, but it is generally understood 
that a great part of the energy dissipation must occur at scales 
that are smaller than the structures that can be resolved by ob- 
servations 100 km), perhaps as small as 10 to 100 m (the 
Kolmogorov turbulent cascade dissipation scale). One of the 
most successful approaches to fill this four-order-of-magnitude 
wide gap is to assume that the statistics obtained at observable 
scales are still valid at smallest scales. The properties of the 
global system, from observable to non-observable scales, can 
then be investigated. This, for example, is the idea underlying 
Hudson's (1991) critical power-law slope of -2 for the distri- 
bution of flare energies. 



) . ) . Statistics from previous observational studies 

Distributions of events. The measurement of the power-law 
slope for the lowest energy flares ha s indeed been a major goal 
of coronal physics in the last decade. IXschwanden et alJ (l2000l) 
has summarized the distributions of event energies that were 
obtained at wavelengths from X rays to ultraviolet (UV), and 
for event energies covering a range of eight orders of magni- 
tude from 10^^ J C'nan oflares") to 10^^ J C ' flares" ). It seems — 
and it is a statement of lAschwanden et all (I200Q|) — that these 



distributions can be matched together to form a unique power- 
law distribution of slope ^ -1.8. 

Spectra. Fourier spectra of the radiance field s in UV and X 
rays have also been obtained (e.g. M artens & Gomez 1992 ; 
iGomez et al]ll993t iBenz et al.lll997MBerghmans et al 
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from Yohkoh/SXT, NIXT and BIT). Most of them display 
power laws (with indices of -1.4 to -4), indicating the proba- 
ble presence of turbulence. 

Structure functions. U sing transverse structure functions, 
Abramen ko et al.l(l2002h have found that the photospheric mag- 
netic field in active regio ns is intermittent (and interm ittency is 
higher for high activitv). [Patsourakos & Viall (l2002h have also 
found intermittency in the time series of radiance (light curves) 
of the lines Ne viii 77.04 nm and O iv 79.02 nm. Moreover, in- 
termittency has been observed in-situ in the solar wind velocity 
and m agnetic field fluctuations LBurlagjll99ll: iMarsch & Tul 
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1.2. Aim of tliis study and selection of the data 

Almost all the statistics of UV observations cited above were 
done on radiance images obtained by pass-band filter imaging 
instruments. Here we want to obtain statistics of the coronal 
turbulence and subsequent heating, and this requires to have 
access to more direct signatures of the basic fields of MHD, 
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namely the 3D magnetic and velocity fields. Some statistics 
of the photospheric m agnetic field have been obtained (cf. 
lAbramenko et al.ll2002h . but they are limited to the line-of- 
sight component. As for coronal magnetic fields, no measure- 
ment has been obtained over a whole solar area. Concerning 
the coronal velocity field, to our knowledge no statistical study 
has yet been done in the framework of turbulence, although ve- 
locity statistics are very important in the study of turbulence. 
To address this subj ect, we need to us e spectroscopic data, and 
SoHO/SUMER (^Wilhelm etalJh995h is still one of the solar 
UV spectrographs which has the best performances currently 
available. Furthermore, a spectrograph also allows us to get 
"purer" data in radiance, because the line radiance is obtained 
by line fitting which removes the continuum background. 

In order to do statistics, we need large amounts of data, 
and also large fields with a good spatial resolution if we were 
to compute Fourier spectra or structure functions over a wide 
range of scales (or frequencies). The resolution of SUMER 

V) allows us to reach scales which are almost the small- 
est scales which are observable by current UV instruments. On 
the other side of the range of scales, scales of the order of the 
size of the Sun can be reached by SUMER when rastering the 
whole image of the Sun, as has been done in the "full-Sun" 
observation program which is detailed in Sect. 12. II 

We also choose to restrict ourselves to quiet regions of the 
Sun. We are looking for statistical results on the turbulent state 
of the corona and on the resulting events of energy dissipa- 
tion that heat the whole corona, active and quiet, and there is a 
priori no reason to exclude specific features like the supergran- 
ulation pattern or active regions. On the other hand, the sta- 
tistical weight of one active region may be too high compared 
to the total amount of data we have, and the characteristics of 
a particular active region may influence the overall statistics. 
For this reason, we exclude data containing active regions from 
our analysis (but we do not perform any particular treatment 
concerning the supergranulation); an analogous study of active 
regions should be the subject of a subsequent work. 

2. The data 

2.1. Observation program and data set 

We use data from a series of 36 observations of the full Sun, 
done in 1996 by one of us (PL). These data, which are listed 
in Table [I] were obt ained by the SUMER spectrograph aboard 
the SoHO satellite (Iwilhelm etal.lh995h . and taken from the 
SoHO archive at MEDOC^ At this time, most of the Sun was 
quiet, and active regions were rare. In each observation of this 
program, 8 rasters in the east- west direction were done with slit 
number 2 (!'' x 30C), and they were accumulated in the north- 
south direction so as to obtain full-Sun images. The exposure 
time was 3 s and the slit moved (in most cases) by 4 elementary 
steps (1.52'0 in the east- west direction. 

Spectra were obtained with detector A of SUMER, but con- 
trary to other observations with this instrument, these spectra 

^ The MEDOC archive is in public access from 
|http://www.medoc-ias.u-psud. fr, 



(or parts of spectra) were not sent to the ground: only 5 param- 
eters of 3 spectral lines were transmitted for each position on 
the Sun (i.e. position of the slit and pixel along the slit). These 
parameters were computed on board by SUMER. During this 
process, some information was of course lost, but this allowed 
us, on the other hand, to get spectroscopic maps of the whole 
Sun at a high resolution (V x 1.5''), with low telemetry use 
(8 • 10^ bits), and quickly (within 9 h). 

The following parameters (or "moments") of spectral line 
profiles were obtained^: 

- (1) peak spectral radiance, (2) Doppler shift, and (3) 
width of the line Svi 93.3 nm (transition 2p^3s^Si/2 - 
2p^3p^P3/2 at 93.340 nm, of mean peak of spectral radiance 
0.57Wm-2sr-inm-i), 

- (4) line radiance (integrated spectral radiance) of the line 
Ly£: (transition ls^Si/2 - 6p^P3/2 at 93.780 nm, of mean 
peak of spectral radiance 1.07 Wm"^sr"%m"^), 

- (5) line radiance of the line Svi 94.4 nm (transition 
2p^3s^Si/2 - 2p^3p^Pi/2 at 94.455 nm, of mean peak of 
spectral radiance 0.29 Wm"^sr"%m"^) 

Hereafter the moments will sometimes be referred to by their 
numbers. 

The Ly6: line is emitted at a formation temperature of 
around 10000 K in the high chromosphere, and the S vi lines 
are emitted around 200 000 K in the transition region between 
the chromosphere and the corona. The width ("moment" 2) is 
computed (on board) as the width of the interval where the line 
profile is higher than half its peak values, and the Doppler shift 
(moment 1) is the position of the center of this interval. The 
data are then compressed aboard SoHO (to one byte per mo- 
ment per pixel), and decompressed on the ground before being 
analyzed (see details of compression in AppendixlAll. 

In addition, some full spectra were recorded so as to check 
for the positions of the spectral windows (used for the on board 
computations) and for the global reliability of the computa- 
tions. One of these is shown in Fig.[l] No specific problem was 
detected. The design of SUMER, in particular the absence of 
cosmic rays impacts on the detector, makes actually the compu- 
tation of line parameters aboard quite reliable, notwithstanding 
the noise problem (see below in Sect. 12.31) . 

2.2. Correction of the data 

As all the data do not come with line profiles, the usual pro- 
cedures of dark current removal and flatfield correction cannot 
be applied. Thus we have to assume that the variations along 
the slit of averages of data for diff'erent dates and positions 
result from systematic instrumental eff'ects and not from so- 
lar structures (actually, large structures like active regions were 
excluded from the computation of the averages). These aver- 
ages, shown in Fig. |2l are then subtracted from the raw data 
(for Doppler shift) or divide the raw data (for radiances), so as 
to get corrected data. For convenience, the units of the resulting 
data will hereafter be named "data units" and denoted du. 

^ Spectroscopic data a nd mean peak of spectral radiances (given for 
quiet Sun) are taken from lCurdTet aL 120011) . 
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Name 


Date 


Start 


End 


Notes 


0401 


1996/04/01 


12:30:43 


17:20:57 


1 


0407 


1996/04/07 


20:24:54 


04:56:18 




0414 


1996/04/14 


01:02:49 


09:34:15 




0418 


1996/04/18 


20:11:39 


04:43:05 




0424 


1996/04/24 


11:08:45 


19:40:12 




0429 


1996/04/29 


01:26:38 


09:58:03 




0504 


1996/05/04 


07:43:09 


16:14:36 




0508 


1996/05/08 


05:09:03 


13:40:28 




0512 


1996/05/12 


23:02:03 


07:33:28 


AR 


0517 


1996/05/17 


00:47:02 


08:32:43 


PI 


0524 


1996/05/24 


08:57:22 


17:28:48 


PI 


0528 


1996/05/28 


12:39:04 


21:10:29 


PI 


0603 


1996/06/03 


21:13:56 


05:45:21 


PI 


0608 


1996/06/08 


21:10:31 


05:41:56 


P2 


0612 


1996/06/12 


16:12:47 


00:44:13 


PI 


0616 


1996/06/16 


22:41:43 


07:13:10 




0621 


1996/06/21 


22:48:43 


07:20:10 


P3 


0701 


1996/07/01 


19:42:26 


04:13:48 


P3 


0706 


1996/07/06 


22:09:36 


06:41:02 


AR P2 


0711 


1996/07/11 


17:44:09 


02:15:35 




0716 


1996/07/16 


19:13:27 


03:44:49 


P2 


0721 


1996/07/21 


18:38:19 


03:09:48 




0726 


1996/07/26 


00:19:04 


08:50:28 




0801 


1996/08/01 


21:07:56 


05:39:21 


AR 


0806 


1996/08/06 


19:14:26 


00:04:04 


1 


0811 


1996/08/11 


16:19:29 


00:50:53 


AR 


0814 


1996/08/14 


06:07:09 


14:05:05 


P3 


0816 


1996/08/16 


20:28:22 


04:59:45 


AR P2 


0821 


1996/08/21 


14:41:18 


23:12:42 


P2 


0828 


1996/08/28 


09:52:13 


12:53:48 


AR P3 


0902 


1996/09/02 


00:11:00 


08:42:22 


P3 


0906 


1996/09/06 


04:23:15 


12:54:41 




0913 


1996/09/13 


21:36:32 


06:07:58 


P3 


0924 


1996/09/24 


15:22:43 


23:53:26 


AR 


0930 


1996/09/30 


22:02:21 


06:33:47 




1005 


1996/10/05 


04:42:48 


13:14:14 





Table 1: List of the observations of the full Sun done in 1996. 
The 4-digit name corresponds to the date of the start of the ob- 
servation (which may end on the following day). Notes: AR, 
active region near the center of the field of view; PI, some 
pointing problems, or missing data, outside of the center of the 
field of view; P2, these problems may aff'ect the center; P3, im- 
portant problems of pointing or missing data; 1 , the raster step 
in the east- west direction is 1.14'' (3 elementary steps) instead 
of 1.52'' (4 elementary steps). 

2.3. Data values and noise 

The context spectra (Fig. [ij, which were taken with an expo- 
sure time of 300 s, show most line profiles with a few hundreds 
of total counts at peak, whereas the full- Sun images were taken 
with a 3 s exposure time. The counts on the detector during 
the full-Sun observations are therefore low, which is confirmed 
by the distributions of the values of radiance moments (1), (4) 
and (5) shown in Fig.|3d), and noise is an actual issue. Monte 
Carlo simulations were performed to estimate the influence of 
counting statistics on the data: we started from Gaussian pro- 
files with background, whose amplitudes were chosen accord- 
ing to the profiles of the context spectra (0.1 counts • pixel"^ 




400 600 800 

Spectral detector pixel 



Figure 1: Top: context spectrum (full detector A) taken on 4 
May 1996 at 7 h 32 min UT. Bottom: total counts along slit be- 
tween pixels 50 and 299. The S vi lines and the Lyman series 
of hydrogen (from Ly 4 to Ly 12) can be seen. 



for the background, and up to 10 counts • pixel" ^ for the line 
profiles), and their results (noise as a function of the line am- 
plitude, for each kind of moment) are shown in Fig. [3l The 
detector noise, of the order of 10""^ counts • s"^ • pixel" \ is not 
taken into account. Taking an amplitude of 4 counts • pixel" ^ as 
an example, the noise (standard deviation of the line parameters 
computed from the Monte Carlo simulations) is: 

- 1.28 counts • pixel"^ for the peak spectral radiance (as mo- 
ment 1), 

- 1 .64 pixels for the line shift (as moment 2), 

- 1.77 pixels for the line width (as moment 3), 

- 8.4 counts for the line radiance (as moments 4 and 5). 

These values, extracted from the plots of Fig.|3l are consistent 
with the sensitivity predicted during the mission preparation 
Jwilhelm et al. 1989 ). 

Among the radiance and line shift data we use, the most 
acute noise problem comes from moment 2 (Doppler shift of 
Svi 93.3 nm): as can be seen in Fig. |4l most of the data are 
contained in the error bars coming from a l-cr evaluation of 
the noise from the Monte-Carlo simulations, thus noise can ac- 
count for most of the dispersion of the data in moment (2). 
Furthermore this moment is saturated due to the limited width 
of the spectral window used to compute the line shift and due to 
the compression algorithms (Appendix |A|), and this saturation 
is still visible in the wings of the distribution of the corrected 
moment (as seen in Fig. |5jd.2)); however, this concerns only 
0.5 % of the pixels and the effect on the statistics is weak. 
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(1) S VI 93.3 nm peak spectral radiance (counts) 




(2) S VI 93.3 nm line shift (pixels) 




(3) S VI 93.3 nm line width (pixels) 



1.65 ^ 
1.60 ^ 
1.55^ 




(4) Ly g radiance (counts) 




(5) S VI 94.4 nm intensity 




Figure 2: Average profiles along the slit over all images (ex- 
cluding specific structures), for each of the 5 moments, in 
counts or pixels. In particular, the well-known odd-even pattern 
of the detector is clearly visible. Note that these cuts represent 
averages around disk center only. 

Peak spectral radiance Shift 




Amplitude (counts) 
Radiance 



40 - 
20 - 



nplitude (counts) 



Amplitude (counts 



Figure 3: Noise (lower curves) and signal (upper curves, when 
relevant) in each pixel as a function of the amplitude of the 
Gaussian profile used for the Monte Carlo simulations, for peak 
spectral radiance, line shift, line width, and line radiance. 

2.4. Construction of full-Sun images 

The 8 rasters for each observation are put together, first by us- 
ing the pointing coordinates from the headers of the data files, 
then by a further adjustment, done so that the limbs on neigh- 
boring rasters fit exactly. Figure ISja) shows the resulting im- 
ages for observation 0721. However, the boundaries between 
rasters can still be seen, mainly because the observation times 




PDF 



2 4 6 8 10 

S VI peak spectral radicnae (counts) 

Figure 4: Joint probability distribution function (PDF) of 
moments 1 (Svi 93.3 nm peak spectral radiance) and 2 
(S VI 93.3 nm line shift), in data units (counts and pixels). The 
horizontal line superimposed on this plot is the median of mo- 
ment (2), and the other curves are the median plus or minus the 
standard deviation of the noise computed by the Monte Carlo 
simulations (Fig.|3ll, as a function of moment (1), assuming that 
the noise in moment (1) is low. 



of adjacent pixels in diff'erent scans can be separated by up 
to 2 hours, and during this time the Sun has rotated by 1.1 ° 
in longitude, or 19'' at disk center as seen by SoHO. For this 
reason, when the statistics imply spatial information, we work 
with data coming from only one raster at a time. This fact does 
not reduce the interest of using these full- Sun images to do 
statistics, as we still have data with a large number of pixels 
produced in a short time. 

2.5. Profiles of the spectral images along the equator 

We produce averaged profiles of the corrected moments along 
the equator of the Sun, in a 200'' -wide band on either side (note 
that here and hereafter we approximate the Sun by a sphere and 
the projection of the solar equator with the solar disk diameter 
which is parallel to the X axis). To do this, we need to (1) take 
into account the variation of the apparent size of the Sun (due 
to the eccentricity of SoHO's orbit around the Sun), which can 
be seen in Fig. |6l the solar X coordinates of the profiles are 
normalized to a common solar diameter before computing the 
average of the profiles; (2) take into account the curvature of the 
limb along the width of the band used to compute the average 
profiles. Thanks to this method, the profiles, shown in Fig.|3b), 
present interesting features, like strong limb brightenings in the 
S VI lines (moments 1 and 5), which can also be seen in the 
images of Fig|3a). The average profiles of these radiances are 
consistent with the theoretical 1 /ji limb brightening of optically 
thin lines (with ji = cos a and a the angle between the normal 
to the solar surface and the line of sight). The fit for moment 1 
(Svi 93.3 nm spectral radiance) gives (0.14 ± 0.05) -h (0.75 ± 



(1) Peak spectral radiance 
S VI 93.3 nm 



(2) Line shift 



(3) Line width 
Svi 93.3 nm 



(4) Radiance 
Ly 



(5) Radiance 
S VI 94.4 nm 




Figure 5: (a) Reconstructed images of the full Sun, on 21 July 1996, for the 5 moments computed by SUMER aboard SoHO. (b) Average cuts along the equator of the full-Sun 
images, for all five moments, (c) Detail of the central field on 21 July 1996, for all five moments, with an aspect ratio of 1. (d) Histograms of the values of the moments on a 
large central field on 21 July 1996. The range of the linear grayscale for each moment is shown on these histograms and is chosen so that 99 % of the pixels fall into this range. 
Note that all the data shown here are the moments corrected according to Sect. 12.21 
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1020 -1000 -98 
Solar X (arcs( 



980 1000 1020 
>olar X (arcsec) 



Figure 6: Profiles of S vi 93.3 nm radiance in a 100-pixel- wide 
band around the equator, as a function of solar X, for all ob- 
servation dates. The units of radiance are data units, and each 
profile is shifted by 100 units from the previous one. 



6 r 





Solar X (arcsec) 

Figure?: Average profile of moment (1) (Svi 93.3 nm peak 
spectral radiance, in data units) along the equator, fitted to the 
l/fi theoretical function. Pixels at the limb and outside of the 
limb, where the optically thin plan-parallel approximation is 
not valid anymore, have been excluded from the plot 



0.03)/yU with;^^ = 5.7 • 10"^ and is shown in Fig.El A limb 
brightening in the Ly s radiance (moment 4), although weaker 
because this line is optically thicker than the S vi lines, has also 
been measured and can be seen in Fig.|3a-b). 

The Svi 93.3 nm line shift (moment 2), also shown in 
Fig. 13b), presents a curvature that can be interpreted as fol- 



lows. Let 6 and (f be respectively the heliographic latitude and 
longitude, with the origin on the line of sight towards the ob- 
server. Assuming that the Sun is a sphere of radius R, rotating at 
an angular velocity fl(6) around a north-south axis perpendic- 
ular to the line of sight^, and assuming that matter flows down 
at a radial velocity u(6, (f), the observed velocity, projected on 
the line of sight, with positive values for downflows, is: 

^los ~ ^(0) ' R cos 6 ' simp -\- u(6, if) ' cos 6 ' cos if (1) 
= Q(0) • X -h u(0, if) • ^1 - (jc2 -h y^)IR^ 

where x = R cos 6 - sin (f and y = RsinO are the heliocentric 
coordinates. The fit of the profile of moment (2) by this func- 
tion (Fig.Ell gives Q = (5±5)' 10"^ du/'' and w = 1.8 ± 0.1 du 
with;^^ =1.9 ' 10"^, where du is the data unit for moment (2). 
Furthermore, a null velocity corresponds to -1.6 ± 0.1 du. 

On the other hand, the unit of line shift given by the de- 
compression routines is one pixel on the detector in the direc- 
tion of spectral dispersion, i.e. 0.0044 nm at a wavelength of 
A = 94 nm, or a 14 kms"^ downflow. We would thus expect 
that Q = 15 • 10"^ du/''. There is a disagreement with the value 
of Q obtained previously, indicating that the fit of solar rotation 
from average equatorial velocity profiles is perhaps not a good 
way to determine the velocity calibration in the case of these 
observations. The origin of this problem may be the low signal 
(see section l23l . which makes it difficult to see the effects of 
solar rotation (expected to be ±2 km • s"^ on the limb). We re- 
tain the value 14 km s"^ for the velocity unit, which seems to be 
the most reliable and which gives physical values for u which 
are more consistent with previous observations in the same line 
dRrekke et alll QQTtlChae et alll 998HPeter & Tudgell 99qh. 

Finally, the Svi 93.3 nm line width (moment 3) is also 
larger near the limb. This increase is easily explained by an 
increase at the limb (due to projection) of the optical depth of 
this optically thin line: the velocity field of the matter emitting 
in this line is averaged over a longer line-of- sight near the limb 
than at disk center. 



3. Statistical analysis 

3.1. Statistics of structures 

Bright structures, which are visible in images, can be identified 
with energy dissipation events. Here we use the same method 
as m lAletti etalJ (E00nL to define structures: one structure is 
a set of pixels whose radiance is above a threshold Lthr, i-e. a 
maximal connex part (cluster) of the set of image pixels with 
radiance above the threshold. Such structures, found using mo- 
ment (1) (Svi 93.3 nm spectral radiance) on 21 July 1996 and 
with a threshold Lthr set to the average radiance plus one stan- 
dard deviation, are outlined on Fig.|9l(top); 2313 of them were 



^ We consider that the heliographic latitude Bq of the disk center 
is 0, i.e. that SoHO is in the equatorial plane; this is an approxima- 
tion due to the angle between this plane and the ecliptic (7.23°), and 
because of the excursions of SoHO's orbit outside the ecliptic. But it 
would only add a second-order contribution to the restriction of Eq. (ij 
to the equator (^ = 0). 
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-500 500 

Solar X (arcsec) 



Figure 8: Average profile along the equator of moment (2) (line 
shift of S VI 93.3 nm, in data units), adjusted to the theoretical 
profile given by Eq.[I] 



found in a central area of 700 x 270 pixels. Since the time dif- 
ference between successi ve images is too long we cannot use 
the time information as in lParnell & Juppl(l2000l) . 

Let Pk = {pj,k, j ^ Jk) be the set of pixels in a structure 
labeled k (Jk is the set of the indices j of these pixels); the radi- 
ance (e.g. moment 1) of the pixel pj^k is called Lj^k- The inten- 
sity 4 of the event of energy dissipation corresponding to the 
structure labeled k is then defined as the integrated radiance on 
the structure, i.e. 4 = YjjeJk ^j,k- For the central field of 21 July 
1996, we obtain a power-law distribution"^ of index - 1 .53+0.03 
over more than two decades, as can be seen on the bottom of 
Fig. |9l This analysis is repeated on the 25 best images, and 
we obtain 40920 structures, whose intensity is distributed as a 
power law of index -1.57+0.01 (Fig.fTol. Compared to the dis- 
tribution for 21 July 1996, this distribution is extended towards 
high intensities because of rare, but intense, events appearing at 
some other dates. The results for the other radiance moments 
are similar and are summarized in Table |2| The error bars are 
the uncertainties on the results of the linear fit of the histograms 
in logarithmic scale, taking into acount Poissonian statistics in 
the counts of events in each histogram bar. However, they don't 
take into acount the uncertainties coming from the selection 
of the fitting range, which could be evaluated to 0.1 approxi- 
mately. 

These slopes a re steeper th a n the slope -1.2 found with a 
similar method by lAletti et al.l ilOOd on SoHO/EIT images. 
On the other hand, they are less stee p than the slopes of the dis- 
tri butions of event energies found bv lAschwanden et all (l2000l) 
or IParnell & JuppI (|2000) for example, who use a diff'erent 
method: they use time information (like time clustering to find 
events), and they use some hypotheses to derive event energies 
from their radiance data. This derivation is highly non-trivial. 



^ The probability distribution functions (PDF) are obtained by 
building an histogram and then dividing the height of the bars by the 
total number of elements and by the width of the current bar. 



Table 2: Slopes of power-law distributions of structures found 
in the maps of radiance moments (Svi 93.3 nm spectral radi- 
ance, Lys radiance and S vi 94.4 nm radiance). Note that the 
given uncertainties do not take into account those coming from 
the choice of the fitting range (see text). 



Obs. 


Moment (1) 


Moment (4) 


Moment (5) 


time(s) 


Svi 93.3 nm 


Lys 


S VI 94.4 nm 


On 0721 


-1.53 + 0.03 


-1.46 + 0.03 


-1.67 + 0.03 


On average 


-1.57 + 0.01 


-1.44 + 0.01 


-1.74 + 0.01 



50- 
0. 



1 Sigma 




ID 100 
Structure rodicinfifi 



1000 



Figure 9: Top: structures (black clusters) determined by a spec- 
tral radiance threshold of S vi 93.3 nm (average plus one stan- 
dard deviation), on a part of the central field, on 21 July 1996. 
Bottom: probability distribution function (PDF) of the spectral 
radiance integrated over these structures, and its fit to a power 
law. The dotted line represents one event per variable-width 
histogram bar and is shown as an indication of the statistical 
noise coming from the construction of the histogram. 



as it depends on the relation between temperature and ther- 
mal energy, on the relation between the area and thickness of 
the emitting structures, on temperature and density diagnostics 
from observations, on the filling factor, etc. Here we have cho- 
sen to remain close to the observable variables, so as to get re- 
sults independent from the assumptions needed to derive event 
energies. We have to keep in mind though that integrated ra- 
diance is not the event energy, and as it is not computed over 
the whole spectrum, it is not even the radiative loss during the 
event. 
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Figure 10: Distribution of the S vi 93.3 nm radiance integrated 
over structures defined as in Fig.|9land over the 25 best images, 
and its fit to a power law. The dotted line represents one event 
per variable- width histogram bar. 



3.2. Fourier spectra 

The omnidirectional Fourier spectrum P(k) of a field is ob- 
tained by computing the integral of the 2D Fourier power spec- 
trum P2D(kx, ky) of the field (which takes into account the dif- 
ferent resolutions in the x and y directions, and using a cosine 
apodization function in both directions) over a ring of radius k: 




Jo 



P(k)dk = kdk de P2D(kcos6,ksme) 



(2) 



In practice, this is done by computing the average of P2D over 
the grid points in 2D Fourier space that fall in the range [k - 
Sk/2, k -h Sk/2] and by multiplying the result by Ink. 



Frequency (orcsec ^) 

Figure 11: Spatial Fourier spectra of moments 1 (S vi 93.3 nm 
spectral radiance) and 2 (S vi 93.3 nm line shift), on a central 
field of the Sun, on 21 July 1996. 

3.3. Structure functions 

For a scalar field a(x), we define the increments S{a(x) = a(x -h 
i) - a(x). The structure function S^(£) of index q is the ^-th 
moment of the increments S^a of the field at scale i.e. : 



S%i)^{\d,a{xr),, 



S VI 93.3 nm spectral radiance. FigureFTTll shows the spatial 
Fourier spectrum of the field of moment 1 (S vi 93.3 nm spec- 
tral radiance) as a function of the spatial frequency / = kjln. 
This spectrum is fitted to a power law of index -1.81+ 0.03 for 

scales lower than 25 Mm (corresponding to a 0.03''"^ spatial S^q^^{£) = S^{£)l {s^{£)\ 
frequency), where a break occurs. This scale is at the order of 
the scale of the supergranulation, and it is the same scale as the 
break in t he Yohkoh/SXT X -ray radiance Fourier spectrum ob- 
tained although the slopes of both parts 
of the spectrum are diff'erent there than in the case of our ob- 
servations. 



where x and x + i go through all the field at a given 
normalized structure functions are: 

ql2 



(3) 



II. The 



(4) 



S VI 93.3 nm line shift. The power-law tail of the spectrum of 
the radiance field indicates that small scales exist, and they 
may be produced by turbulence. However, the direct measure- 
ment of the spectrum of the line-of- sight velocity (or line shift 
of S VI 93.3 nm) — a physical parameter easier to compare to 
spectra predicted by theories of turbulence — does not yield 
satisfying results: these spectra are flat, as can be seen from 
Fig.[n]2 for 21 July 1996, implying that the noise in the maps 
of this parameter (moment 2) cannot be ignored. 



In particular, S^^^^ii) is called the flatnes s, F{t). 

In phenomenological turbulence like [Kolmogoro^ (Il94lh 
or Iroshnikov-Kraichnan (llroshnikovll 19631: iKraichnanI 
the structure functions of the velocity are power laws of index 
with oc q, and this leads to normalized structure functions 
that are constant as a function of the scale £. When higher-order 
structure functions, in particular the flatness, depend on it is 
a deviation from the above-mentioned classical phenomenolo- 
gies of turbulence and the field is intermittent. 

Svi 93.3 nm spectral radiance. The normalized structure 
functions for the scalar field of moment 1 (S vi 93.3 nm spec- 
tral radiance) averaged over 20 observations of the quiet Sun at 
disk center are shown in Fig. [T2| They clearly rise when the 
spatial scales gets smaller, especially when they get smaller 
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Figure 12: Normalized structure functions of index 1 to 5 (from 
bottom to top) for the field of moment (1) (S vi 93.3 nm spectral 
radiance), averaged over 20 observations of the quiet Sun at 
disk center; the thick line is the flatness. 

than 15 Mm. This rise is characteristic of correlations in the 
field and of intermittency, and can also be seen in the normal- 
ized structure functions of the other radiance moments (of Ly s 
and S VI 94.4 nm). However, when studying turbulent fields, the 
fields of which the intermittency is of primary interest are the 
magnetic and velocity fields. Consequently we now examine 
the latter. 

S VI 93.3 nm line shift. When the field is a vector field a(x), it 
is usually projected longitudinally when computing the struc- 
ture functions (in Eq. 6ea(x) is replaced by Sea(x) • 
In o ur case, as well a s in magnetic field observation like 
lAbr amenko et al] (l2Q02l) . we have no choice but to do this pro- 
jection along the line-of-sight unitary vector ei^^: 

SH{) = {\(v(x + e)-v(x)).eiJl^^^^^^__^ (5) 

The resulting average normalized structure functions are shown 
in top of Fig. El Because of the noise (which is not spatially 
correlated) in the fields of moment (2) and perhaps because of 
other eff'ects discussed later, these structure functions are much 
flatter than in the case of moment (1), even if a small rise at the 
small scales can still be seen. However, when plotting the flat- 
nesses of the field of moment (2) at all dates used to compute 
the average (Fig. El bottom), one sees that the rise of the flat- 
ness is systematic, and this gives a higher level of confidence in 
the rise observed on average. It seems thus that we can observe 
intermittency in the velocity field of the S vi 93.3 nm line. 

4. Discussion 

Many difficulties arise when attempting to compute the statis- 
tics presented here, mainly because of trade-offs between the 
need to get large fields at high spatial resolution during a lim- 
ited time (during which the Sun has not changed much), and 
to: 



Flatness, moment (2) 



10 




10 100 
6x (arcsec) 



Figure 13: Top: normalized structure functions of index 1 to 5 
(bottom to top) for the field of moment (2) (S vi 93.3 nm line 
shift), averaged over 16 observations of the quiet Sun in the 
disk center; the thick line is the ffatness. Bottom: flatness of the 
field of moment (2) at the dates used to compute the average 
(which is shown with the thick line). 

- send all the required data to the ground with a given band- 
width (or to keep them temporarily in the limited SUMER 
memory) 

- have long enough exposure times so as to have enough sig- 
nal. 

Because of the first trade-off we need to compute the moments 
of the lines on board and compress them by limiting their dy- 
namics (see Appendix |A||, and the second trade-off' is at the 
origin of the noise problem, which is particularly important for 
the line shift data (moment 2). But with the current instruments 
there is no other alternative: we think that the data set we used 
here is one of the most adapted to the purpose. 

Despite the difficulties, the tentative results presented here 
support the presence in the transition region of small scales 
created by turbulence and the intermittent nature of turbulent 
fields. To our knowledge this is the first determination of the 
intermittency of the velocity field in an ultraviolet line. These 
small scales, at which dissipation is more efficient, may play a 
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role in the heating of the plasma of the solar upper atmosphere, 
particularly of the corona (if confirmed by similar observations 
in hotter lines). 

Improved observations should be obtained with more sen- 
sitive spectrographs, such as the ones from Solar B and Solar 
Orbiter, which will allow us to make progress by increas- 
ing spatial resolution while also increasing counting statistics 
and thus allow, in particular, a better determination of the ve- 
locity field. However, there will still be compromises to be 
made on the path to "ideal" observations which would have 
less noise, better spatial and temporal resolution, and which 
would include lines emitted at hotter tempratures further up 
in the corona, closer to where the heating is believed to take 
place. One step further could be attained by using novel ul- 
traviolet imagin g spectrograph designs such as the MOSES 
spectro-imager (iKankelborg & Thomasll200l|) or an Irn aging 
UV Fourier Transform Spectrometer JMillard et al]l2004 . 

Appendix A: Some details about the compression 
algorithm 

The need for sending to the ground the data of several full- Sun 
images (one per moment) at a high resolution, in a limited time, 
and with a limited bandwidth (10.5 kbps), forced us not only to 
compute the moments on board ("compression" method num- 
ber 17 of the SUMER Operations Guide) but also to limit the 
number of bits used for one moment in one pixel. The principle 
of this second "compression" is simply to reduce the dynamics 
to a 8 -bit unsigned byte and to saturate high (and low, in the 
case of line shifts) values. 

To perform this, the function near is defined by the code 
corresponding to: 

near(jc) = byte (min (x -h 0.5, 255)) (A.l) 

where min (•, •) returns the smallest of its both arguments and 
byte (•) transforms its argument into a one-byte integer lower 
approximate value. 

If Ml , . . . , M5 are the 5 moments computed by the com- 
pression scheme number 17 for one pixel of the image, then 
the values transmitted to the ground are: 





= near (255 Ml /Ci) 


(A.2) 


M2 


= near (max (10 M2 + 128,0)) 


(A.3) 




= near (10 M3) 


(A.4) 


m; 


= near (255 M4/C2) 


(A.5) 


^5 


= near (255 M5/C3) 


(A.6) 



where Ci is the COMPAR/ parameter given in the headers of the 
FITS file. 

The inverse operation is performed on the ground to recover 
approximations of the original values of Mi , . . . , M5 . 

The on-board computation of the moments and the further 
compression of the data were necessary to get the results pre- 
sented in this paper from SoHO/SUMER data. However, fu- 
ture observations of this kind should learn from the limitations 
found here: even if they could still benefit from such a method, 
there is a need for more signal in the original spectral data, and 



for a compression scheme that allows a better dynamics (like a 
lossless compression scheme applied on data with more than 8 
bits). 
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